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ABSTRACT 

The universality and linearity of the Leavitt law are hypotheses commonly adopted 
in studies of galaxy distances using Cephcid variables as standard candles. In order 
to test these hypotheses, we obtain slopes of the Leavitt law using linear regressions 
of fundamental-mode Cepheids observed by the Optical Gravitational Lensing Exper- 
iment project in the Magellanic Clouds. We find that the slopes in FT-bands and in 
the Wcscnhcit index behave exponentially, indicating non-linearity. We also find that 
the slopes obtained using long-period Cepheids can be considered as universal in the 
y/-bands, but not in the Wesenheit index. 
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1 INTRODUCTION 

During the first decade of the twentieth century, ob- 
servations of variable stars from the Small Magellanic 
Cloud (SMC) galaxy taken with the 24-inch Bruce tele- 
scope at Harvard Boyden Station in Arequipa, Peru, 
conducted to Henrietta Leavitt to make an important 
discovery: the correlation between the logarithm of the 
pulsation periods of Cephei d variables and their magnitudes 
l|Leavitt fc Pickering! Il912t l. In honour of its discoverer, 
this correlation, commonly called Period-Luminosity (PL) 
relatio n, has been renamed rece ntly as the Leavitt law 
(LL) (|Freedman fc Madord lioiol ). Originally, the size of 
the sample used by Leavitt to establish the PL relation 
was statistically small: only 25 Cepheids. Thirteen years 
later, the number of SMC Cepheids used to establish 
the PL relation increased by more than a factor of four 
(|Shaplev. Yamamoto fc Wilson! Il925l ). More recently, mi- 
crolensing projects, as the Massive Compact Halo Objects 
(MACHO) and the Optical Gravitational Lensing Experi- 
ment (OGLE), discovered hundreds of Cepheids, as a result 
of monitoring systematically the bars of the Magellanic 
Clouds. Based on these observations, it was possible to 
determine, with high precision, periods and mean VI 
magnitudes of Cepheids, and the slope and zero-point 
values o f the LL for t he La rge Magellanic Cloud (LMC) 
galaxy (jlldalski et all Il999ah . These advances confirmed 
the LL as an useful modern tool to determine extragalactic 
distances, ti ed to an assumed LM C distance modulus of 
18.50 mag (|Freedman et al.l l200lf ). as shown by several 
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studies in the last years, but in particular through the 
work of the Araucaria Project by measuring distances 
to selec ted galaxies of the Local Group and the Scu lptor 
Group (|Pietrzvnski et al.|[200l . IPietrzyfiski et al.ll20l61 . and 
references therein). 

In order to measure distances at galaxies through the 
LL, two important hypotheses have to be adopted: its 
universality and linearity. Universality means that the 
slope of the LL is independent (or weakly dependent) on 
metallicity, implying that this value in a photometric band 
is the same for all galaxies. Linearity means that there 
is a linear correlation between mean magnitudes and the 
logarithm of the pulsation periods of Cepheids. 
The influence of metallicity on the PL relation has been dis- 
cussed by sev e ral gr oups. Studying some sel e cted g alaxies, 
lUdalski et ail (|200lt ) and IPietrzyfiski et al.l (|2004 ) found 
strong evidence in favour of the universal slope of the PL 
relation in optical bands, between the m etallicity range of 
-1.0 d ex (I C1613) up to -0.3 dex (LMC) . IPietrzvnski et ail 
|2004l ) and iBono et al] \201(t ) showed that the use of the 
Wesenheit index, Wi = I — 1.55(V — J), minimises not only 
the effects of reddening bu t also those of metallicity on the 
PL relations. iGieren et all (120051 ) and iFouque et all (|2007l ) 
showed that the slopes of the PL relations in V7-bands and 
in the Wi index do not change significantly in the Milky 
Way and LMC galaxies. 

Although the universality and linearity of the PL relation 
have been assumed in distance determination studies, there 
are works showing that these hypotheses are not always 
truthf ul. With respect to the non-universality, IStorm et al.l 
(|201lD , using the infrared surface brightness technique, 
found that the slopes and zero points of the PL relations of 
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the LMC and Milky Way present a significant dependence 
on metallicity in the Wesenheit index, a weak dependence 
on metallicity in the optical V7-bands, and do not present 
any dependence on metallicity in K-haxiA. 
With respect to non-linearity, the Russian astronomer 
Boris Kukarkin found in the thirties that the PL relations 
of external galaxies presen t a bre ak around 10 days, as it 
was mentioned by Ferriiel d 19691) a nd references therein. 
iTammann. Sandage fc Reindll (|2003h found a break of the 
Period-Colour (PC) relat ion at 10 days in the LMC galaxy. 
iKanbur fe Ngeowl |2004) confirmed these break based on 
the statistical _F-test. 

From a study of the LMC Cepheids, IKanbur et"afl (|2007l l 
detected non-linearity of the LL. This was found when using 
the testimator method, that determines if there is a slope 
change in a nth subset with respect to a smoothed slope 
obtained in all previous subsets. The smoothing process is 
applied in order to avoid the outliers influence on the slopes 
computed. Applying this method, those authors detected 
a significant change in the slope between subsets of period 
around 10 days. 

iKoen. Kanbur fc Ngeowl (|2007h presented a study of 
non-linearity of the PL relation and the Period-Luminosity- 
Colour (PLC) relation using two non-parametric meth- 
ods: linear regression residuals and additive models. 
These authors found again that the PL relation is non- 
linear and presents a break around 10 days. Possible 
causes for the non-linearit y of th e LL were suggested by 
IKanbur. Ngeow fc Buchlerl and IKanbur fc Ngeowl 

(l200rt). 

iNgeow et all (|2012h studied the PL relation at individual 
phases of the pulsation cycle of Cepheids (the multi-phase 
approach) and found again a break around 10 days. They 
showed that the behaviour of the LL at mean light is the 
average of the behaviour at all phases. This result reinforce 
the use of the PL relation at mean light in the distance 
determination context. Readers interested in more details 
of these statistical methods are referred to the papers cited 
above. 

Despite the results of these studies, universality and 
linearity of the LL remain to be adopted. The slope 
and zero point of the LL in optical V7-bands have 
been established with a high degree of accura cy from 
funda mental- mode Cepheids in the LMC galaxy bv lUdalskil 
When the ordinary least-squares (OLS) regression 
was applied to OGLE-II data of Cepheids, over a range of 
logP from 0.4 up to 1.5 for the LMC, and from 0.4 up to 1.7 
for the SMC, Udalski realised that the standard deviation 
of the residuals of the LMC data was almost two times 
smaller than that of the SMC. Despite that the number of 
Cepheids in the SMC doubles that of the LMC, the greater 
dispersion of the SMC PL relation is caused by the spatial 
distribution of Cepheids in the galaxy bar, whose thickness 
is placed alo ng the line-of-sight with a typical depth of 
~ 0.25 mag (|Harris fc Zaritskvl I2006T ). This fact, and the 
manner in which the PL relation of the LMC populates for 
periods greater than 2.5 days , led to adopt the LMC slope 
value as a universal quantity ([Udalski 2 000T ). 

These facts and results we have mentioned above, motivate 
some questions behind this paper: what is the slope of the 



LL from a larger sample of Cepheids? More specifically, 
what is the slope if the Cepheids in the LMC and SMC are 
considered as a single sample? 

In order to test the universality and linearity hypotheses 
of the LL and answer the previous questions, this paper 
is structured as follows: a description of OGLE-II Cepheid 
data of the Magellanic Clouds is presented in the second 
section. The extension of the OLS regression, that we 
develop to calculate the slope of the LL from Cepheids in 
the Magellanic Clouds simultaneously, is found in the third 
section. Tests of universality and linearity of the LL, as also 
our main results and a discussion of them, are presented in 
the fourth section. Finally, our main conclusions are given 
in the fifth section. 



2 DATA 

In a study of the PL relation in optical bands using 
LMC OGLE-III fundamental-mode Cepheids, INgeow et al.l 
( 2009) found a discrepancy in the slopes computed by them 
and those reported by other authors, due to the different 
number of Cepheids used in each study. For this reason, 
to avoid effects of the sample size, in our study we use 
OGLE-II data to establish a d irect compariso n between our 
slopes and those obtained by lUda lski (2000), that remain 
as the accepted universal values. 

The OGLE II photometric data for the Magel- 
lanic Clouds were obtained betwe e n 1997 and 1999 , 
(|Udalski. Kubiak fc Szvmahskil 1 19971 : ISzvmariskil I2005T ). 
High quality BVI observations of hundreds of variable 
stars along the galaxies bars were collected with the 1.3-m 
Warsaw telescope, a t Las Campanas Observatory, Chile 
(|Udalski et al.lll999al ). A number of epochs between 120 to 
360 was obtained in the /-band and between 15 to 40 in 
the _BF-bands. A search for periodicity in the photometric 
time-series of th e stars was made, using the a nalysis of vari- 
ance algorithm (|Schwarzenberg-Czernv|[l989f ). The Cepheid 
variables were selected based on visual inspection of their 
light cur ves and their locati on on the Colour-Magnitude 
diagram (|Udalski et aljl999h l . To select fundamental-mode 
Cepheids, an analysis of Fourier coefficients was performed. 
From the OGLE-II on-line archive we select a number of 
745 and 1287 fundamental-mode Cepheids, belonging to 
the LMC and SMC, respectively. Their dereddened mean 
magnitude s were computed onc e the extinction effect was 
measured (|Udalski et alj Il999bl , lUdalski et ail Il999ch . To 
compute the Ai extinction coefficient, it was determined 
the /-band magnitude of Red Clump stars in many lines- 
of-sight. Differences of the observed /-band magnitude of 
Red Clump stars were assumed as differences in the Ai ex- 
tinction coefficient. To compute the colour excess E(B — V) 
and the extinction coefficient Ay in each line-of-sight, they 
adopt ed the reddening law of ISchlegel. Finkbeiner fc Davisl 
(|l998ft . 



3 LINEAR REGRESSIONS 

To test the universality and linearity hypotheses of the LL 
we are interested in obtaining the slopes of the PL relations 
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in optical bands without adopting the values of the slopes 
reported bv lUdalskil (|200Ch . With this purpose, we apply two 
different approaches of linear fits: the OLS regression and the 
multiple least-squares (MLS) regression, that is an extension 
of the OLS regression developed by us. The MLS regression 
is a similar approach to the testimator method described in 
the introduction section, but with these method we want 
to find the slope of the PL relation when the mathemati- 
cal union of data sets of Cepheids from different galaxies is 
used, under the hypotheses that all data sets share the same 
slope (universality hypothesis), but each one has a different 
zero point. 

Before presenting the MLS regression, we briefly recall how 
the OLS regression works, emphasising on some useful re- 
sults. 



3.1 Ordinary Least-Squares Regression 

llsobe et al.1 l|l990t ) stated that the OLS regression should be 
used to predict the value of one variable from the measure- 
ment of another, under the following conditions: 

(i) The nature of the linear regression scatter is not under- 
stood, and this scatter is always greater than the errors of 
the measured variables. 

(ii) It is well established that one variable is the cause and 
the other is the effect. 

As a first approximation, this is the case of the LL: we want 
to predict the mean magnitudes of Cepheids from their mea- 
sured periods. To do this, it is necessary to use the OLS 
regression over a set of N Cepheids, each one of them char- 
acterised by two variables: logPfc and (m) k , where the k 
sub-index takes values between 1 to N. The logarithm of 
the Cepheid period is log P^ , and in the context of this 
study, {m) k is the dereddened mean magnitude in a pho- 
tometric band. The residuals of the measurements of dered- 
dened mean magnitudes, d k , are defined as the square of the 
differences between the observed values of these mean mag- 
nitudes and the values predicted by the OLS regression, as 
follows: 

d k = {{m) h -r 1 \ogP k -0' 2 - (1) 
The standard deviation of the residuals is given by: 



T {rn) 



(N-2) : 



(2) 



where S is the sum of the residuals dk. The expressions of the 
zero point, £, and the slope, r\, can be found by minimising 
dk- The well know results are: 
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(3) 
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It is possible to express the standard deviations of the slope 
and zero point in term s of the standard deviation of the 
residuals (|Tavlorlll982T ). as follows: 
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where A is given by the following equation: 



A = iV^(logP fe ) 2 - E lo S P ' 



(6) 



3.2 Multiple Least-Squares Regression 

Let G be a set of galaxies, each one having a different num- 
ber Ni of Cepheids. We want to find an analytical expression 
for the slope of the PL relation resulting from the mathe- 
matical union of Cepheid data of all G galaxies, under the 
assumption that all data sets share the same slope r\. We 
call this value the common slope. We also want to find the 
analytical expression for the zero point corresponding to 
each galaxy, under the assumption that each galaxy is at a 
different distance. For the Ith galaxy, the PL relation can be 
described as: 

{m) kl =ri\ogPu+ii, k = 1,2,..., Ni I = 1,2,..., G. (8) 

(m) k , is the predicted dereddened mean magnitude of the 
kth Cepheid belonging to the Ith galaxy, from its measured 
period P k i . In analogy with equation fl} , we define the resid- 
uals of the measurements of dereddened mean magnitudes 
for the Ith galaxy, as: 



Ni 



A = E (( m )ki -v^ogPm - 6) s 



(9) 



We also define the sum of the residuals for all galaxies as: 

G 

T(^)=E^- ( 10 ) 

!=1 

To find the values of the common slope r\ and the zero points 
£i, T must be minimised with respect to each one of the 1+G 
variables (77, £(): 

on & 

By solving these equations, the expression for the zero points 
£1 of each galaxy is obtained as: 

&4 EHr-iEWiL (12) 

; \fe=l fc=l / 

which is a similar expression to the equation Q, given by 
the OLS regression, but with 77 being the common slope 
which is independent of the zero points and that is given 
by the following equation: 



n = 



EHi Efcir (log^) 2 - Ef =1 i(E£r logPfcO 2 
where 



;=i \fe=i / \fc=i 



. (13) 



(14) 



The equations ([3]) and Q can be obtained of the equations 
(j 121) and (|13|l considering a unitary set of galaxies. In the 
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case of the MLS regression it is possible to obtain an expres- 
sion for the standard deviation of the residuals generalising 
the equation ([2]), taking into account that it is necessary to 
include the data points of all PL relations. This expression 
for the standard deviation of the residual magnitudes can 
be written as: 



r <™> 



T 



(15) 



The standard deviations of the common slope rj and the zero 
points can be obtained by generalising the equations $5$ 
and ©: 



(16) 



(17) 



where A is given by equation 0. 

The above expressions are useful in order to test the univer- 
sality and linearity hypotheses of the PL relation. In order 
to achieve that, we apply the MLS regression on Cepheids of 
the LMC and SMC, as a statistically sample of Cepheids. In 
this case, the I index runs up to 2. We will refer hereafter to 
the data set formed by the mathematical union of the LMC 
Cepheids with the SMC Cepheids as the LMC + SMC set. 
At this point, it is important summarise the main charac- 
teristics of the linear regression methods exposed. The OLS 
regression allows to obtain the slope and zero point of the 
PL relation of a single galaxy. By applying this method to 
the SMC, an increasing of the dispersion of the PL relation 
caused by the distribution of the Cepheids along the galaxy 
bar is noted. In order to obtain simultaneously the slope and 
zero points of the PL relations using the LMC and SMC 
Cepheids, and damping the effect of the dispersion caused 
by the SMC Cepheids, we use the MLS regression over the 
LMC + SMC set. 

This method allows us to test the universality of the LL 
in VJ-bands making a study different to the traditional 
ones: we want to establish if the slopes obtained from the 
LMC + SMC set are consistent with the slopes reported by 
Udalski (2000) for the LMC galaxy, despite of the metallicity 
differences of the Magellanic Clouds. 

3.3 Testing the Multiple Least-Squares 
Regression 

Before applying the MLS regression to the observational 
data, we perform some tests of the MLS regression over sim- 
ulated data sets in order to find the common slope and the 
zero points of each set and verify the validity of the results of 
this regression method. First, ten sets of PL relations, each 
one with 10000 data points, are randomly generated. Each 
data set has a different zero point and the same slope. Then 
the MLS regression is applied, obtaining the same values 
of slope and zero points used to generate the PL relations. 
Then, random noise is added to the generated periods, and 
the MLS regression is applied again. The slopes and zero 
points obtained using the MLS regression are the ones ex- 
pected. In these tests we do not observe trends in the residu- 
als of the mean magnitudes with respect to log P. Finally, we 
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Figure 1. Percentage difference of removed stars vs. 
threshold value. Results for the V /-bands and the Wj index are 
shown in the upper, middle and bottom panels, respectively. D 
presents a local minimum in a = 2.5. Triangles (black) correspond 
to the LMC Cepheid data, squares (red) show the SMC Cepheid 
data and circles (blue) correspond to the LMC + SMC set. A 
colour-version of this figure is available in the on-line journal. 



increase the number of sampling points by one order of mag- 
nitude. As expected, we observe a corresponding decrease of 
the standard deviations of slopes and zero points, as well as 
in the dispersion of the residuals of the mean magnitudes. 
As the MLS regression works properly with synthetic data, 
we apply it to the observational data sets of Cepheids in the 
Magellanic Clouds. 



3.4 Removing outliers in the PL Relations 

Before to obtain the slopes and zero points of the LL us- 
ing the OLS and MLS regressions, the outliers must be dis- 
carded. To remove them, we find the standard deviation of 
the residual magnitudes, <Ji m ) and <Ti m ) , given by equations 
(01 and (|15|l . Then, points with a dispersion larger than 
CKX/ m ) and c«T/ m \ , are removed. In order to find the opti- 
mum threshold a, we study the behaviour of the percentage 
difference of removed stars vs. a, selecting a between 2 and 
3, with steps of 0.1. Denoting by RS(a) the percentage of 
removed stars with threshold a, and RS(a — 0.1) the per- 
centage of removed stars with threshold a — 0.1, we define 
the percentage difference as: D = RS(a) — RS(a — 0.1). 
Fig. [TJ shows a plot of D vs. a for the LMC, SMC, and the 
LMC + SMC set, in the V7-bands and the Wi index. By 
examining this figure, it is clear that the local minimum in 
the studied range lies in a — 2.5. For this minimum value 
there are equal numbers of rejected points between the ad- 
jacent threshold values of 2.4 and 2.5, indicating that it is 
the optimum threshold value. It is important to note that 
this threshold value is the same for t he three bands studied , 
and is equal to the value reported bv lUdalski et all (|l999bh . 
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Figure 2. Histograms of log P for Magellanic Clouds 
Cepheids observed by the OGLE-II project. The LMC 

sample presents a sudden decrease of detected Cepheids for 
logP < 0.4 (upper panel). For the SMC and the LMC + SMC 
set, the number of detected Cepheids decreases for logP < 0.2 
(middle and bottom panel). 



3.5 Selecting Appropriate Ranges of Periods 

In order to obtain slopes and zero points of the PL relations, 
it is necessary to establish the range of periods that should 
be selected to apply the linear regressions. Since Cepheids 
with shortest periods are faintest, they may have been diffi- 
cult to detect with the 1.3-m Warsaw telescope, used by the 
OGLE project. Cepheids with longest periods, that are too 
bright, were over-exposed. Fig. [5] shows a strong decrease in 
the number of detected Cepheids in the LMC galaxy about 
logP = 0.4, and in the SMC galaxy and the LMC + SMC 
set about logP = 0.2. The saturation level of the CCD used 
by the OGLE-II project, determined the upper limit of de- 
tectio n in logP = 1.5 for t he LMC and logP = 1.7 for the 
SMC l|Udalski et al.lll999al ). 



3.6 Discarding Cepheids with logP > 1.5 

We also study Cepheids with periods larger than the mea- 
sured by the OGLE-II project (logP > 1.5), using photo- 
metric data from the third phase of the All Sky Automated 
Survey (ASAS-3, |Poimariski|[2012T ). This data base contains 
photometry from almost ten years of continuous observa- 
tions. The ASAS-3 VI magnitudes are near to the Johnson's 
standard system, but they were computed without include 
colour terms in their calibration^ The number of Cepheids 
that we find with periods between 30 and 100 days is 19 
and 7 for the LMC and SMC, respectively. Our search for 
pe riodicity is made using the an alysis of variance algorithm 
by ISchwarzenberg-Czernvl l|l989t ). 

According to the photometric accuracy and the phase cov- 
erage, we try fit each light curve by Fourier series of orders 
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ranging from 3 to 6. The mean magnitudes ar e comp uted 
following the method described bv lGieren et all (|2004T ). For 
each one of the ASAS-3 Cepheids, the best fit was obtained 
using a Fourier series of 4i/i order. The mean magnitudes 
obtained from other orders, produce very small effects on 
the LL. We determine that the changes induced on the zero 
points and slopes by using these mean magnitudes are of the 
order of hundredths. 

When the ASAS-3 Cepheids are included in the PL relations, 
they are placed too far from the trend of Cepheids observed 
by the OGLE-II project. There are two main causes that 
could explain this behaviour. The first one is the difficulty to 
correct by extinction, because these Cepheids are placed far 
from the galaxies bars, where the extinction has been deter- 
mined with less accuracy. The second one is the large errors 
(up to several tens of magnitudes) found in the transforma- 
tion of the instr umental magnitudes to the stand ard system 
for some stars (|Poimanski fc Macieiews ki 2005). Based on 
these causes, we decide not to include the ASAS-3 Cepheids 
in our study. 

iKarczmarek et all (|201ll 'l calculated the slope of the PL 
relation from 65 LMC Cepheids, selected from the ASAS 
catalogue of variable stars. They found a slope value of 
Vasas = —2.366 ± 0.166 mag/dex for the V-band. This 
value di f fers in more than 10a from the value reported by 
lUdalskil (|2000l ) and from the value of the slope calculated 
by us with MLS regression (see Table Q]). This fact supports 
our decision to exclude these Cepheids of the sample. 



4 RESULTS AND DISCUSSION 

In order to test the universality and linearity hypotheses 
of the LL, we calculate the slopes of the PL relations for 
the Magellanic Clouds Cepheids through two different ap- 
proaches. By using the OLS regression it is obtained the 
slope of the PL relation of each galaxy. By using the MLS 
regression we obtain the slope value of the LL from Cepheids 
belonging to the LMC + SMC set. 



4.1 Testing the Universality Hypothesis of the 
Leavitt Law 

The universality hypothesis of the LL can be tested by show- 
ing that the value of the slope in a photometric band is the 
same regardless the metallicity of the sample of Cepheids 
used to derive it. We compute the slope of the Leavitt law 
for the LMC + SMC set using the MLS regression, assum- 
ing the linearity of this law in the Vl-bsaids and in the Wi 
index, and using the ranges of periods mentioned in sub- 
section 3.5, 0.2 < logP < 1.7 for the SMC galaxy and the 
LMC+SMC set, and 0.4 < log P < 1.5 for the LMC galaxy. 
Table[T]presents the slopes and their st andard deviatio ns ob- 
tained by us (MLS regression) and by lUdalskil (|2000T I (OLS 
regression ). It is worthy of mentioning that the results ob- 
tained bv iudalskil (|2000h have been exactly re-obtained in 
this work. In order to establish a comparison between our 
results and the reported by Udalski, the slope values t)mls 
are computed using dereddened VI magnitudes. The first 
column in Table [T] gives the photometric band in which the 
LL is studied, the second and third columns give the com- 
mon slopes and their standard deviations obtained by using 
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Table 1. Slopes obtained testing the universality hypoth- 
esi s. Columns 2 a nd 4 give the slope values obtained by us and 
by lUdalskil j200Cf ). respectively. A detailed explanation of each 
column is given in the text. 



Band 


VMLS 


a n 


VOLS 


a n 


Nlmc 


NsmC 


V 


-2.806 


0.028 


-2.775 


0.031 


700 


882 


i 


-3.024 


0.020 


-2.977 


0.021 


705 


857 


Wj 


-3.364 


0.012 


-3.300 


0.011 


729 


821 



Table 2. Asymptotes fi obtained testing the linearity hy- 
pothesis. A complete explanation of each column is given in the 
text. 





LMC 


SMC 


LMC+SMC 


Band 


H 


A/i 




A/i 




V 


-2.726 


2.751 


-2.759 


0.025 


-2.767 0.062 


1 


-2.909 


2.130 


-3.008 


0.029 


-2.993 0.040 


Wi 


-3.311 


0.005 


-3.376 


0.013 


-3.358 0.009 



the equations (|13[) and (|16[) . Columns fourth and fifth give 
the slopes and the i r stan dard deviations for the LMC ob- 
tained by lUdalskl l|200(j ). applying the equations Q and 
((51 . Columns sixth and seventh give the total number of 
Cepheids in the LMC and SMC galaxies that we use with 
the MLS regression, after rejecting outliers. Our computed 
common slope t/mls of the LL is consistent at a level of 
1.0a and 2. 2a in V- and J -bands, respectively, with those 
reported by lUdalskil (|200d ). This implies that the slope of 
the PL relation in the U/-bands is the same independently of 
the metallicity of the Magellanic Clouds. On the other hand, 
th e slope Vails i n the Wi index is not consistent with that 
of lUdalskil | |200Ct ). These facts indicate that the LL is uni- 
versal in the W-bands but not in the Wi index. The slope 
obtained in J-band is not as accurate as the reported by 
Udalski, probably due to blending effects reinforced by the 
bar geometry of the SMC galaxy. If Cepheids are blended 
with red stars, this may change their mean / magnitudes 
increasing the dispersion of the PL relation. 

4.2 Testing the Linearity Hypothesis of the 
Leavitt Law 

The linearity hypothesis of the LL can be tested by show- 
ing that the value of the slope in a photometric band is 
the same regardless of the range of periods selected to de- 
rive it. In order to compute the slope, we select the range 
of periods between a lower limit and an upper limit. The 
lower limit is useful to avoid the effect of Malmquist bias 
on the slope of the LL, that makes the slope value too shal- 
low. This effect is present when the faintest Cepheids in a 
galax y are very close to the c utoff magnitude of the photom- 
etry (|Pietrzvnski et aj]|201Ch . According to the discussion 
in subsection 3.5, we take the lower limit for the range of 
periods at logP = 0.4 for the LMC, and logP = 0.2 for 
the SMC and the LMC + SMC set. We select the upper 
period limit of the Cepheids ranging from logP = 0.5 up to 
log P = 1.5, with steps of 0.1, for the LMC, SMC and the 
LMC + SMC set. 
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Figure 3. Slope vs. logarithm of the upper limit for range 
of periods. Triangles and stars represent the slopes obtained 
through the OLS regression for the LMC and SMC, respectively. 
Filled circles show the common slope obtained with the MLS 
regression, for the LMC + SMC set. Filled triangles are identi- 
fying the LMC outliers data, see text for explanation. Dashed 
lines (black for the LMC, red for the SMC and blue for the 
LMC + SMC set) are the exponential fits given by equation 
l|18|l. Upper, middle and bottom panels show the behaviour of 
the slope of the LL for the U7-bands and for the Wi index, re- 
spectively. In each panel are show the slope errors obtained by 
using the MLS regression over the LMC + SMC set. Large error 
bars corresponds to slopes computed with few data between the 
range of periods. The slope errors corresponding to the LMC and 
SMC, are similar to the errors of the LMC + SMC set. They are 
omitted to appreciate clearly the figure. A colour-version of this 
figure is available in the on-line journal. 



Fig. [3] shows the slope values for the PL relations as a func- 
tion of the logarithm of the upper limit of period. Each one 
of these slope values for LMC and SMC samples was ob- 
tained by using OLS regression, and the slope values for the 
LMC + SMC set were obtained using the MLS regression. 
The behaviour of the obtained slopes is well fitted by an 
exponential function (dashed lines in Fig. [3}: 

T) = Ae -/31ogP +Ati (lg) 

where the A, /3, and fj, parameters are obtained using a non- 
linear fitting with the IRAF task nfitld. The A parameter 
gives information about the stretching or squeezing of the 
function and the reflection about the horizontal axis. The /3 
parameter gives the reflection about the vertical axis. The 
asymptote n can be interpreted as follows: if the Cepheid 
population of a galaxy is observed, including all Cepheids 
with periods up to 100 days, and their mean magnitudes 
and periods are used in order to obtain the slope of the PL 
relation, this slope value approaches to (i. 
Inspecting the Fig. [3] it is clear that the LMC LL presents a 
break around 10 days, in U/-bands. This result is in agree- 
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ment w ith those reported by Tammann. Sandage fc Reindll 
J2003h. iKanbur fc Ngeowl h2004 ), iKanbur et al.l (|2007h and 
iKoen. Kanbur fc Ngeowl (|2007l ). 

The LMC slopes in 1/7-bands behave anomalously for short- 
periods (logP < 1.0, black triangles in Fig. [3}, in the sense 
that they move away from the exponential function which 
adjusts the data trend for long-periods (logP > 1.0). For 
this reason short-period points for the LMC are considered 
to be outliers to fit the data by the equation (|18p . On the 
other hand, the obtained slopes by the OLS regression for 
the LMC are shallower than the slopes of the SMC for short- 
periods. The causes of this behaviour are still under inves- 
tigation by us. 

The behaviour of the slopes shown in Fig. [3] provides clear 
evidence of the non-linearity of the LL in the V7-bands and 
in the Wi index for the Cepheids belonging to the Magel- 
lanic Clouds. This result about the non-linearity of the LL 
is in concordance with the c orresponding non- l ineari ty of 
the PC relation, reported by Kanbur fc Ngeowl (|2004l ) and 
IKoen. Kanbur fc Ngeowl (|2007T ) for the LMC Cepheids. Be- 
sides, due to that the PC and the PL relations are projec- 
tions of the fundamental PLC relation (|Madore fc Freedmanl 
1991). our result of non-universality of the PL relation sug- 
gests that the PC relation should be non-universal in the Wi 
index. In Table [2] we present the values of the asymptote ji 
and their errors obtained by using a non-linear fit with the 
IRAF task nfitld. The errors of the asymptote ix, A/x, are 
computed taking into account the slope errors. The first col- 
umn shows the photometric band where the LL is studied. 
The second and third columns give for the LMC the values 
of n and A/x. The large Afj, values for the LMC reflect that, 
in that case, the 11 values are obtained using about one half 
of the points that are used in other fits. The following pairs 
of columns in Table [2] give the same information that the 
first two columns, but for the SMC and the LMC + SMC 
set. 

Despite the non-linearity of the PL relation and its break 
around 10 days (for the LMC), its use to determine distance 
is not significantly impacted by these effects due to the fol- 
lowing reason: 

The calculation of distance modulus of a galaxy, relative 
to the LMC, is performed from the difference between zero 
points of the PL relations of the galaxy and the LMC. These 
zero points are function of the slope of the PL relation as 
shown the equations and (j 12 j) . 

The slope of the PL relation is influenced strongly by long- 
period Cepheids. In a galaxy these Cepheids are less nu- 
merous than short-period Cepheids and are distributed in a 
wide range of periods. Therefore, a small number of them 
can produce a significant change in the slope value more 
than a large number of short-period Cepheids. 
The effects of non-linearity of the PL relation can be esti- 
mated computing the deviation of the slope from the value 
of xx. Inspecting Fig. [21 it is seen that the bigger is the period 
range used to derive the slope of the PL relation, the smaller 
is the effect of non-linearity. For example, in V-band the de- 
viation of the slope obtained with MLS regression from the 
value of /x, in the period range from 2.5 to 10 days, is 3a, 

while for a period range from 2.5 to 25 days , is 1.1a . 

The slopes of the PL relation calculated by lUdalskH (|2000l ) 
were obtained with Cepheids up to 32 days (LMC) and 50 
days (SMC), for which the effect of non-linearity obtained 



in the present work is in the order of one sigma in V-band. 
Based on these facts our results reinforce the use of the LL 
to determine distances, independently of the non-linearity 
and the break of the PL relation. 

It can be seen from Table [2] that the values of the asymp- 
totes fj, for the LMC + SMC set, in VI -bunds, agree in 
less than 1.0<r with the value s for the slopes of the PL rela- 
tions found bv lUdalskil (|2000l ) (see Table □ fourth and fifth 
columns). The /x values obtained for the LM C and SMC ar e 
consistent with the values of r\ published bv lUdalskil ||2000D . 
up to 1.6cr for the V-band, and 3.2a for the /-band. Also, 
it can be seen from Table [2] a statistically significant differ- 
ence between the (x values of the LMC and SMC in the Wi 
index, and also betwe en the /x value of each galaxy with the 
r\ value published by lUda lski (20Q3). These results suggest 
that the slope /x can be considered as universal in V7-bands 
but not in the Wi index. In a forthcoming paper we will 
be applying the OLS and MLS regressions in galaxies other 
than the Magellanic Clouds, in order to test the universality 
and linearity of the LL in these metallicity environments. 



5 CONCLUSIONS 

In this work we test the universality and linearity hypothe- 
ses of the LL, using the sample of Cepheids belonging to the 
Magellanic Clouds, observed in 1/7-bands by the OGLE-II 
project. In order to develop these tests, we compute the 
slope values of the LL using two different approaches. One 
of them makes a mathematical union of the Cepheid data 
of the LMC and SMC galaxies to find a common slope ap- 
plying the MLS regression. The other one obtains the slope 
of the PL relation applying the OLS regression on a single 
galaxy: the LMC or the SMC. 

The test of the universality hypothesis lead us to obtain 
a common slope for the LMC + SMC set using the MLS 
regression. Ou r values are consistent with the reported by 
lUdalskil (|2000l ) for the LMC at a level of l.Ocr and 2.2cr in 
V— and /—bands, respectively, and inconsistent in the Wi 
index. Our results suggest a strong dependence on metallic- 
ity of the slope for the Wi index, and a weak dependence 
on metallicity in the o ptical V7-bands, in agreement with 
the results reported bv lStorm et al.l |201ll ). reaffirming the 
universality of the LL in the W-bands and showing again 
its non-universality in the Wi index. 

The test of the linearity hypothesis lead us to find that 
the values of the slopes of the LL behave exponentially as 
the range of periods increases, giving clear evidence of non- 
linearity. In particular, we find that the LMC LL presents 
a break around 10 days, according to t h e results reported 

DW 



a break around 10 days, accordin g to t h e results re portc 
by Tammann. Sandage fc Reindll (2003), Kanb ur fc Ngeo 



Kanbur et al.l ( 20071 ) and Koen. K anbur & Ngeow 



(2004). 

( 2007). A clear explanation of the behaviour of the LMC 
PL relation, in the range of short-periods is unknown un- 
til now. Our result about the non-linearity of the LL is 
in concordance with the corre sponding non-linearity o f the 
LMC PC relation, reported bylKanbur fc Ngeowl (|2004l ) and 
IKoen. Kanbur fc Ngeowl (120071 ) , and suggests that in the Wi 
index the PC relation should be non-universal. 
Despite of the non-linearity of the LL, we find that the 
asymptote xx has value s in V /-bands t hat are consistent with 
the slopes reported bv lUdalskil (|2000l ) for the LMC, but in- 
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consistent in the Wi index; therefore, in V/-bands, fj, can 
be cons i dered as universal and the LMC slopes given by 
lUdalskil l|200Cl ) remain appropriate to measure extragalac- 
tic distances. However, it is necessary to take into account 
that the universality of the LL has been tested only in a 
few galaxies; therefore it will be necessary to study the PL 
relation in many galaxies other than the Magellanic Clouds, 
in order to test the universality of the LL in different metal- 
licity environments. 
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